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Ll 
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LLL IIIIIIIII BBBBBBBBBBBB RRRRRRRRRRRR TTTTTTTTTTTTTTT «LLL Ll 
LLL HII BBBBBBBBBBBB RRRRRRRRRRRR TTTTTTTTTTTTITTT «LLL + 
LLL LITT BBB RRRRRRRRRRRR TTTTTTTTTITITITT «LLL L 
LLL 111 BBB BBB RRR RRR TTT LLL | Ll 
LLL II] BB BBB RRR RRR TTT LLL LI 
LLL 11] BBB BBB RRR RRR TTT LLL | Ll 
LLL II] BBB BBB RRR RRR TTT LLL | Ll 
LLL Il] BBB RRR RRR TTT LLL Ll 
LLL Il] BBB BBB RRR RRR TTT LLL Ll 
LLL III B RRRRRRRRRRRR TTT LLL Ll 
LLL Ill BBBBBBBBBBBB RRRRRRRRRRRR TTT LLL LI 
LLL Ill 8 RRRRRRRRRRRR TTT LLL ul 
LLL III BRA BBB RRR RRR TTT LLL 
LLL II] BBB RRR RRR TTT LLL Ll 
LLL II] BER BBB RRR RRR TTT LLL Ll 
LLL II] BRB BBB RRR RRR TTT LLL Ll 
LLL III ERB BBB RRR RRR TTT LLL Ll 
LLL Ill BSB BB RRR RRR TTT LLL Ll 
LLELLLLLLLLLLLLL IIIIIIIII BBBBBBBBBBBB RRR RRR TTT LELLLLLLLLLLLLL Ll 
LLELLLLLLLLLLLLL TITIII111 BBBBBBBBBBBB RRR RRR TTT LI 
LLELLLLLLLLLLLLL LIIIIIII BBBBBBBBBBBB RRR RRR TTT LELLLLLLLLLLLLL 
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LL II1111 BBBBBBBB PPPPPPPP 000000 LL YY YY HH HH . 
LL I11111 BBBBBBBB PPPPPPPP 000000 LL YY YY HH He ; 
LL I] BB BB PP PP 00 oO LL YY YY HH HH 3 
LL I] 8B BB PP PP 00 oo LL YY YY HH HH ; 
LL II BB BB PP PP 00 0O LL YY YY HH HH 3 
LL Il BB BB PP PP 00 oO LL yY YY HH HH 3 
LL Il 68888888 PPPPPPPP 00 oo LL YY HHHHAMHH HHH 3 
LL li 8B BB PPPPPPPP 00 0O LL YY HHHHHMM AHH, F 
LL II 8B BB PP 00 oO LL YY HH HM 3 
LL Il BB BB PP 00 00 LL YY Hid HH 3 
LL II BB BB PP 00 oo LL YY Hid HH ceee : 
LL II BB BB PP 00 0o LL YY HH HH eane : 
LLLLLLLLLL IIIIII BBBBBBBB PP 000000 LLLLLLLLLL YY HH HH “ens : 
LLLLLLLLLL III] BBBBBBBB PP 000000 LLLLLLLLLL YY HH HH anaw : 
LL IIIII1 SSSSSSSS 3 
LL HII! SSSSSSSS 3 
LL II ss 3 
LL Il SS 3 
LL II ss 3 
LL II SS é 
LL I] SSSSSS 3 
LL II SSSSSS 3 
LL Il ss 

LL I] ss 

LL I] Ss 

LL II ss 

LLLLLLLLLL III111 SSSSSSSS 

LLLLLLLLLL HII] SSSSSSSS 
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j Cueeeecesecencerceccceseeseseqnnseecenerenensneeceseeeseseseeenseceeeesoeees 
COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 


te Ys A ag UNDER A LICENSE AND MAY BE USED AND oor ace 
N_ ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WiTH THE 
ION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
Pp 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
aORPORATI ou NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR _ RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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; LABABAARSAASAASASLALALALAAASAALAASESESESESEEE LESSEE ES ELESEASSEL ESAS EEE ELE ESS 
+e FACILITY: General Utility Library 
ABSTRACT: 
Perform H floating point polynomial calculation. 
ENVIRONMENT: User Mode, AST Reentrant 


AUTHOR: Steven B. Lionel, CREATION DATE: 31-July-1979 
MODIFIED BY: 
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- Use local handler to insure that ex 


as statuses are resignalled. 
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v04- p 
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H floating polynomial calculat 16-SEP-1984 00:16:51 VAX/VMS M 
NS =" Se Cae te 


.SBTTL DECLARATIONS | 
: INCLUDE FILES: | 
SCHEDEF 
$SSDEF 


co 
m 
oO 
am) 
> 
z— 
ro 
“4 
—-3 
Oo 


Condition handling symbols 
System symbols 


> EXTERNAL DECLARATIONS: 
.EXTRN LIBS$SIG_TO_RET 


Library routine to convert 
a signal to an error return 
to caller of LISS Ve 

RO = signalled condition 


Sete Se Se 
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> MACROS: 


PEGI 


MEAN" OOONOU LS Wwr-00 
. . 


6 
68 ; EQUATED SYMBOLS: 
0 e 
00000004 71 arg = 4 3 argument 
000000 7 degree = | : eegree of polynomial 
000000 7 coeff = 1 ; address of coefficient 
000000 74 result = 16 F 


result of polynomial 
table 


77 : OWN STORAGE: 


-PSECT _LIBSCODE PIC, USR, CON, REL, LCL, SHR, - 
EXE, RD, NOWRT, LONG 


—-9ooO 
one 
SOOOCOOCOCOCOOOOOOOOOOOOOCOCOO 


SCooQoooooocoo 
SooQoooooooo 
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: > PSECT DECLARATIONS: 
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LIBSPOLYH = Pertormef oating polynomial ae 99:09:84 LIBRTL.SRCJLIBPOLYH.MAR; 1 


7 os -SBTTL LIBSPOLYH = Perform floating polynomial 
: FUNCTIONAL DESCRIPTION: 


4 88°01 VY provides the functionality of the VAX hardware 
nstruction POLYH to high Level language users. 


The third operand points to a table (array) of H_floating 
coefficients. The coefficient of 
the highest order term of the polynomial is pointed to 


Ooo 


— 


e 
0 39 : by the table address operand, i.e. the first table element. 
$88 38 3 The table is specified with lower order coefficients stored 
B86 18 : at increasing addresses. 
$008 1 } ; The evaluation is carried out by Horner's method, and the 
0000 1 § ; result is stored at the location pointed to by the fourth 
4 32 ; operand. The result computed is: 
0000 105; if d = degree 
0000 196 ; and x = ar 
0000 107 | result = ChOJ+x*(CC1J+x#(CC2I+ ... x#CCd])) 
0000 109; The unsigned word degree operand specifies the highest 
944 119 3 numbered coefficient to participate in the evaluation. 
0000 116 : For further detail, refer to the VAX-11 Architecture 
4 aE : Handbook for the description of POLYx. 
9000 113 : CALLING SEQUENCE: 
0000 117; status.wic.v = LIBSPOLYH (arg.rh.r, degree.rw.r, coeff.rh.ra, 
44 13 : result.wh.r 
9000 0 : INPUT PARAMETERS: 
$900 1 § : arh.rh.r - argument, ‘'x'' in polynomial 
0000 123; degree.rw.r - degree of pol ynontal (GEQ 0) 
i448 : : ; coeff.rh.ra - table of coefficients, H_floating 
0000 1 6 > IMPLICIT INPUTS: 
000 127; 
at aes 
99 0 ; OUTPUT PARAMETERS: 
09 j : result.wh.r - result of calculation 
000 134 : IMPLICIT OUTPUTS: 
00 135; 
09 ! § : NONE 
8 > FUNCTION VALUE: 
9 140 ; SS$_NORMAL - successful execution 
141 ; SS$_FLTOVF - floating overflow 
146 3 SS$_FLTUND - pees ing under flow 
145 ; SS$_ROPRAND - reserved operand, see VAX Architecture 


Page 
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133 ; manual for more details 
146 ; SIDE EFFECTS: | 
138 : ALL other exceptions are resignalled. 
1 0 i=. | 
134 | 
403c 00 136 ENTRY LIBSPOLYH, “M<IV,R2,R3,R4,R5> ; Entry point, eqeble int. ovf. 
38 137 ; and save R2, R3, R4, RS 
6D 17°AF QE 09 133 MOVAB B*“*HANDLER, (FP) ; Enable local handler to process 
B00e 128 i; exceptions 
oc BC 08 BC 04 BC 75FD 0006 158 POLYH @arg(AP), - ; perform polynomial 
OO0OE 159 adegree(AP), ; trap on exception to 
OOOE 160 acoeff (AP) 3; handler which will 
QOOOE 161 ; unwind a return erver 
OO0E 16¢ ; condition in RO t 
443 197 : caller of LiBsPoL YD. 
10 BC 50 7DFD O03 192 mMOovO RO, @result (AP) 3; return value 
50 01 9A oot 19 MOVZBL #1, RO ; success status code 
04 3018 169 RET 3 return 
oo17 «17h 
0017 176 HANDLER: 
0000 0017 17 -WORD 0 
Be 
0019 176: If the exception is one of the documented exceptions for this routine, 
0019 773s eatt t1Bssl TO_RET to return it as a status. Otherwise, resignal. 
opie 178 ; = Also, “restenti Tf the depth is not zero. 
O19 180 | 
50 08 AC 0 $019 181 MOVL CHFSL_MCHARGLST(AP), RO ; Get mechanism vector address 
08 AO D5 0010 86 TSTL HF $L_MCH_DEPTH(RO) ; Is depth zero? 
41 is 3 » 3 BNEQ 0$ ; If not, resignal 
51 O04 AC OD 2 «(184 MOVL CHF ft “B1GMRGL ST te) R1 ; Get signal vector address 
50 04 Al DO 0026 185 MOVL cares NAME(R1), RO ; Get signalled condition 
048C 8F 50 81 OO2ZA 1 § CMPW R PSS$- FLTOVF ; Compare conditions 
eA 13 0 F 4 Bea. «=s«108 ; If it matches, don't resignal | 
049C 8F 9 " 1 : 3 ne . § #SS$_FLTUND 
0454 BF 0 Bt 8 139 Hay . § #SS$_ROPRAND | 
O4ce BF 9G BI O05F 198 CHP RO, #SS$_FLTUND_F | 
04 Al 049C 4 7 rt, 132 ppv els. ; gS _FLTUND. CHFSL_SIG_NAME(R1) ; Change fault code to trap code | 
0484 8F 0 8} 33 136 5$ eee s § #SS$_FLTOVF _F | 
4 Al O048C BF 3C 198 MOVZ2WL #SSS$_FLTOVF, CHFSL_SIG_NAME(R1) ; 
00000000'GF 6 FA 58 199 10S: CALLG (AP)? G*LIB$SIG_TO-RET ; Return signal as a status 


1-602." 


beret * Hence teben' sh Asch “ECRRSIGRE ALISA UMGIRE SESE ARSCER ess OP hy 


lo 
er 
50 0918 8F Be i Re poyeu. #SS$_RESIGNAL, RO ; Resignal condition 
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LIBSPOLYH 
Symbol table 


HF SL_MCH 
CH SL ~SLGRRGLST 
CH PSL ~S1G_NAME 


C 
DEGREE 
HANDLER 


LIBSPOLYH 
LIB$SIG_TO_RET 
RE SUL 
5S$~ROPRAND 


PSECT name 


Initialization 


The _ working se 
21679 bytes (4 


ARG 
ed oy -MCHARGLST 
DEPTH 


Te processing 
eee | table sort 


Symbol table output 
Psect synopsis output 
Cross-reference output 
Assembler run totals 


Page faults 


e 
- Perform H floating polynomial ssleviot 16-SEP AX/VMS Pa 
need ae ane eo ve 


ro V04-00 
REIL IBPOLYH. MAR; 1 


b+ 4 
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00 RG 


gererenre 


ooo 
oror 


Allocation PSECT No. Attributes 

00000000 ( } oF ?-} NOPIC USR CON ABS _ LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
00000000 0:) . NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
00000069 <« 105.) 5°} PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


were ware ama meno eae. eae + 


cr Tine Elepsed Tine 

13g GBs08 88-07 88:89:01 

18, 8o:00:00:59  00:00:02:08 

54 00: 00:00:37 00:00:08-20 
5 00: 00:00-09 00:09:00.0 
§ 00:00:00.0 90:00:00.0 
385 00:00:04.14 00:00:22.5 


Limit was 1050 pages. 
pages) of virtual memory were used to buffer the intermediate code. 
There were 30 pages of symbol table space aliocas d to hold 42 

source Lines were read in Pass 1, produci ng 
9 pages of virtual memory were used to define 8 macros. 


non-local ia local symbols. 


object records in Pass 
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! \ Macro Library Statistics : 


) v0. 
| 


Macro Library name Macros defined 
~$255SDUAZ8: [SYSLIBJSTARLET.MLB; 2 ates ee 
486 GETS were required to define 5 macros. 
There were no errors, warnings or information messages. 
| 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=L1IS$:LIBPOLYH/OBJ=OBJ$:LIBPOLYH MSRC$:LIBPOLYH/UPDATE=(ENH$:LIBPOLYH) 
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